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Abstract 

Dilepton production in nucleus-nucleus collisions at 1-2 A- GeV as well as in elemen- 
tary pp and pd reactions is studied within the microscopic HSD transport approach 
which includes the off-shell dynamics of vector mesons explicitly. The study ad- 
dresses additionally the production of 7r° and rj mesons since their Dalitz decays 
provide a sizeable contribution to the dilepton invariant mass spectra up to about 
0.5 GeV. Our transport results agree with the TAPS experimental data on tt and rj 
multiplicities in l2 C + 12 C collisions from 0.8 to 2 A GeV. We find that the 'DLS- 
puzzle' - which stands for a theoretical underestimation of the e + e~ yield in the 
mass range from 0.2 to 0.5 GeV in 12 C + 12 C and 40 Ca + 40 Ca collisions - may be 
solved when incorporating a stronger bremsstrahlung contribution in line with re- 
cent OBE calculations. Moreover, the HSD results with 'enhanced' bremsstrahlung 
cross sections agree very well with the HADES experimental data for the dilepton 
mass spectra for 12 C+ 12 C at 1 and 2 A-GeV, especially when including a collisional 
broadening in the vector-meson spectral functions. Detailed predictions for dilep- 
ton spectra from pp and pn/pd reactions at 1.25 GeV, 2.2 GeV and 3.5 GeV are 
presented which will allow to verify /falsify the larger bremsstrahlung contributions 
from the experimental side in the near future. 
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1 Introduction 



The theory of quantum-chromo-dynamics (QCD) describes hadrons as many-body 
bound or resonant states of partonic constituents. While the properties of hadrons are 
rather well known in free space (embedded in a nonperturbative QCD vacuum) the 
mass and lifetimes of hadrons in a baryonic and/or mesonic environment are subject 
of current research in order to achieve a better understanding of the strong interaction 
and the nature of confinement. In this context the modification of hadron properties 
in nuclear matter is of fundamental interest (cf. Refs. [1.2 3 4 5 6 78 9 1Q]1 1 1 II 1 2 j for 
the leading concepts) since QCD sum rules [2][o]fT5] as well as QCD inspired effective 
Lagrangian models [T|4][DXrf51l5|10|llirT2^ predict significant changes e.g. of the vector 
mesons (p, u and 0) with the nuclear density pn and/or temperature T. 



A more direct evidence for the modification of vector mesons has been obtained from 
the enhanced production of lepton pairs above known sources in nucleus-nucleus 
collisions at SPS energies [i"4"][T5] . As proposed by Li, Ko, and Brown [16J and Ko 



et al. [T7| the observed enhancement in the invariant mass range 0.3 < M < 0.7 
GeV might be due to a shift of the p-meson mass following Brown/Rho scaling pQ 
or the Hatsuda and Lee sum rule prediction [2]. The microscopic transport studies 
in Refs. [T8|19|20|21|22j for these systems have given support for this interpretation 
[T6fl7f26|r2~7] . On the other hand also more conventional approaches that describe a 
melting of the p-meson in the medium due to the strong hadronic coupling (along the 
lines of Refs. [3" f 4" | 6 f 8 | 12]) have been found to be compatible with the early CERES 
data [8lfT8ll23T24ll25l. 



This ambiguous situation has been clarified to some extent by the NA60 Collabo- 
ration since the invariant mass spectra for p + p~ pairs from In+In collisions at 158 
A-GeV clearly favored the 'melting p' scenario [28]. Also the latest data from the 
CERES Collaboration (with enhanced mass resolution) [2S] show a preference for the 
'melting p' scenario which is in line with more recent theoretical studies on dilepton 
production at SPS energies [301131 .32 33 34 35 36137] . 

Dileptons have also been measured in heavy-ion collisions at the BEVALAC by the 
DLS Collaboration [38|39] at incident energies that are two orders-of-magnitude lower 
than that at the SPS. The first published spectra at 1 A-GeV [38] (based on a limi- 
ted data set) have been consistent with the results from transport model calculations 
that include pn bremsstrahlung, 7r°, n and A Dalitz decay and pion-pion 



annihilation. However, in 1997 the DLS Collaboration released a new set of data 
based on the full data sample and an improved analysis, which showed a consider- 
able increase in the dilepton yield: more than a factor of five above the early DLS 
data [38] and the corresponding theoretical results [4"0~|41f42|T4'3"] . This discrepancy 
remained in the transport calculations even after including contributions from the 
decay of p and uo mesons that are produced directly from nucleon-nucleon and pion- 
nucleon scattering in the early reaction phase [2"2~|T4'4"] . With an in- medium p spectral 
function, as that used in Ref. [23] for dilepton production from heavy-ion collisions at 
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SPS energies, dileptons from the decay of both directly produced p's and pion-pion 
annihilation have been considered, and a factor of about two enhancement has been 
obtained in the theoretical studies compared to the case of a free p-spectral func- 
tion. Furthermore, in Ref. [45] the alternative scenario of a dropping p-meson mass 
and its influence on the properties of the iV(1520) resonance has been investigated. 
Indeed, an incorporation of such medium effects lead to an enhancement of the p- 
meson yield, however, was not sufficient to explain the DLS data. Since independent 
transport calculations by HSD (BUU) and UrQMD underestimated the DLS data 
for C + C and Ca + Ca at 1A GeV by roughly the same amount these findings have 
led to the denotation 'DLS-puzzle' in 1999 [22|44|I45| which persists now by about a 
decade. 

As has been shown in 2003 by the transport analysis of the Tubingen group [47] 
also alternative scenarios for the in-medium modification, i.e. a possible decoherence 
between the intermediate mesonic states in the vector resonance decay, increases 
the dilepton yield. However, still the region about M ~ 0.3 GeV was underestimated 
whereas the yield in the vicinity of the p-meson peak was overestimated (especially for 
C + C collisons). Thus there is no consistent explanation for the DLS data so far and 
questions came up about the validity of the DLS data or a possible misinterpretation 
of the DLS acceptance cuts. 

In order to address the in-medium modifications of vector-mesons from the experi- 
mental side the HADES spectrometer has been built [5T] which allows to study e + e~ 
pair production in a much wider acceptance region for elementary pp, pd reactions as 
well as it A, pA or even AA collisions up to about 8 A- GeV. During the construction 
phase of HADES there has been a lot of theoretical activity providing various pre- 
dictions for these reactions at SIS energies [21,48,49,52,53 54,55 56 57]. Meanwhile 
the HADES Collaboration has presented first spectra [BD1l61f62|[D r 3] and the question 
is: what do these data tell us? 

The answer to the questions raised above is nontrivial due to the nonequilibrium 
nature of these reactions and transport models have to be incorporated to disentangle 
the various sources that contribute to the final dilepton spectra seen experimentally. 
In this study we aim at contributing to this task employing an up-to-date relativistic 
transport model (HSD) that incorporates the relevant off-shell dynamics of the vector 
mesons. Compared to our earlier studies in Refs. [2~T1l23ll44|r4"o] a couple of extensions 
have been implemented such as 

• off-shell dynamics for vector mesons - according to Refs. [64,65j - and an explicit 
(nonperturbative) propagation of p, u and also at SIS energies 

• extension of the LUND string model to include 'modified' spectral functions for 
the resonances in the string decays 

• novel cross sections for r] production in line with more recent experimental data 

• novel NN bremsstrahlung cross sections according to recent OBE calculations 

• an extended set of spectral functions for the vector mesons p, uj, <fi. 
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We note that the first off-shell dynamical study for vector mesons and dilepton 
decays has been presented in Ref. [48j for pA reactions and been picked up later by 
the Rossendorf group [57] also incorporating the off-shell equations from Refs. [64165] . 
The off-shell dynamics is particularly important for resonances with a rather long 
lifetime in vacuum but strongly decreasing lifetime in the nuclear medium (especially 
uj and (f) mesons). 

The outline of the paper is as follows: After a brief description of the transport 
model and its elementary channels in Section 2 we will discuss tt° and r\ production 
in elementary as well as AA reactions in comparison to available data (Section 3) in 
order to have independent constraints on the leptonic sources. Section 4 is devoted 
to dilepton production in pp and pd reactions in confrontation with spectra from 
the DLS while in Section 5 a comprehensive study of 12 C + 12 C and 40 Ca + 40 Ca 
collisions is performed employing different spectral functions for the vector mesons. In 
addition we will present predictions for future projects of the HADES Collaboration 
like Ar + KCl at 1.75 A- GeV or pp, pd reactions at 1.25 GeV, 2.2 GeV and 3.5 GeV. 
A summary and discussion closes this work in Section 6. 



2 Description of the model 



Our analysis is carried out within the HSD transport model [20|21|66] - based on 
covariant self energies for the baryons [58] - that has been used for the description 
of pA and AA collisions from SIS to RHIC energies. We recall that in the HSD 
approach nucleons, A's, N*(1440), N*(1535), A, £ and £* hyperons, E's, S*'s and 
$Ts as well as their antiparticles are included on the baryonic side whereas the 0~ and 
1~ octet states are incorporated in the mesonic sector. Inelastic baryon-baryon (and 
meson-baryon) collisions with energies above y/s th ~ 2.6 GeV (and \fs th ~ 2.3 GeV) 
are described by the Fritiof string model [67] whereas low energy hadron-hadron 
collisions are modeled in line with experimental cross sections. 

The dilepton production by a (baryonic or mesonic) resonance R decay can be 
schematically presented in the following way: 



BB^RX (1) 

mB -> RX (2) 

R -> e + e~X : (3) 

R — > ml, m — > e + e~X, (4) 

R -> R'X, R' -> e + e-X, (5) 

i.e. in a first step a resonance R might be produced in baryon-baryon (BB) or meson- 
baryon (mB) collisions ([T]), Then this resonance can couple to dileptons directly 
(j2J) (e.g., Dalitz decay of the A resonance: A — > e + e~N) or decays to a meson m 
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(+ baryon) or in (jlj) produce dileptons via direct decays (p, uj) or Dalitz decays 
(7T°, 77, uj). The resonance R might also decay into another resonance R' ((Sj) which 
later produces dileptons via Dalitz decay. Note, that in the combined model the 
final particles - which couple to dileptons - can be produced also via non-resonant 
mechanisms, i.e. 'background' channels at low and intermediate energies or string 
decay at high energies [59] . 

The electromagnetic part of all conventional dilepton sources - n°,n,uj Dalitz de- 
cays, direct decay of vector mesons p, u and <fi - are treated as described in detail 
in Ref. [68J- where dilepton production in pp and pd reactions has been studied. 
Modifications - relative to Ref. [68] - are related to the Dalitz decay of baryonic 
resonances and especially the strength of the pp and pn bremsstrahlung since calcu- 
lations by Kaptari and Kampfer in 2006 (69] indicated that the latter channels might 
have been severely underestimated in our previous studies on dilepton production at 
SIS energies. In detail: 

- For the Dalitz decays of the baryonic resonances we adopt the parametrizations 
from Ernst et al. [22] (eqs. (9) to (13)) which are also incorporated in the PLUTO 
simulation program of the HADES Collaboration |70j . 

- For the bremsstrahlung channels in pp and pn reactions we adjust our previous 
expressions in order to match the recent results from the OBE model calculations in 
Ref. [69] (cf. Section 2.5). 



2. 1 Vector-meson spectral functions 

In order to explore the influence of in-medium effects on the vector-meson spectral 
functions we incorporate the effect of collisional broadening (as in Refs. jl5|72] ). i.e. 
the vector meson width has been implemented as: 

r* v (M, \p\,p N ) = T v (M) + T coll (M, \p\,p N ). (6) 

Here Ty(M) is the total width of the vector mesons (V = p,uj) in the vacuum. For 
the p meson we use 



r,(M)^r_(M) = r [£j f(m) (7) 

(M 2 -Am 2 ) 1 ' 2 (M 2 -Ami) 1 ' 2 
q = , qo = ^ • 

In ([7j) Mo is the vacuum pole mass of the vector meson spectral function, F(M) is a 
formfactor taken from Ref. [8] as 
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Fig. 1. The ratio R(M) = a(e + e~ tt+tt") / 'cr(e+ e~ 
corresponds to the result with the formfactor F(M) ([8]) while the dashed (blue) line is 
without F(M). The compilation of the experimental data is taken from Ref. |73j . 



F(M) 



' 2A 2 + M 2 \ 2 
2A 2 + M 2 J 



(8) 



with a cut-off parameter A = 3.1 GeV. This formfactor was introduced in Ref. [8] in 
order to describe the e + e~ experimental data with better accuracy. In Fig. [1] we show 
the ratio R(M) = a(e + e~ — > n + n^)/cr(e + e~ — > p + p~) which is proportional to the 
pion formfactor F n (M) except for a kinematical factor. The solid line corresponds 
to the result with the formfactor F(M) (JHJ) while the dashed line is without F(M). 
The compilation of the experimental data is taken from Ref. [73] . As seen from Fig. 
[T]the effect of the formfactor F(M) ([8]) is visible only at high M. 



For the ui meson a constant total vacuum width is used: T u 
is a narrow resonance in vacuum. 



r^(Mo), since the uj 



The collisional width in ([6]) is approximated as 



T coll (M, \p\, p N ) = 7 p N < v a% >^ a coll ^. (9) 

Po 

Here v = \p\/E; p, E are the velocity, 3-momentum and energy of the vector meson 
in the rest frame of the nucleon current and 7 2 = 1/(1 — v 2 ). Furthermore, pn is the 
nuclear density and ay f N the meson-nucleon total cross section. 

In order to simplify the actual calculations for dilepton production the coefficient 
a co u has been extracted in the HSD transport calculations from the vector-meson 
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collision rate in C + C and Ca + Ca reactions (at 1 and 2 A- GeV) as a function of the 
density p^. In case of the p meson the collision rate is dominated by the absorption 
channels pN — > tcN or pN — > An — > ttttN at SIS energies. Also the reactions 
P + tt <-> a\ are incorporated but of minor importance here. The numerical results for 
FcoIi(pn) then have been divided by Pn / Po to fix the coefficient a co u in Qj. We obtain 
OLcoii ~ 150 MeV for the p and a co u ~ 70 MeV for u mesons which are consistent with 
Ref . pjl] . In this way the average effects of collisional broadening are incorporated in 
accordance with the transport calculations and allow for an explicit representation 
of the vector-meson spectral functions versus the nuclear density (see below). We 
point out, however, that only the average width of the vector mesons is consistent 
with the dynamical evolution of these particles in the transport calculations and not 
the specific shape of their spectral functions. 

In order to explore the observable consequences of vector meson mass shifts at finite 
nuclear density - as indicated by the CBELSA-TAPS data [75] for the u meson - 
the in-medium vector meson pole masses are modeled (optionally) according to the 
Hatsuda and Lee [2] or Brown/Rho scaling [1] as 



M °*w = jrr^y (10) 

where pN is the nuclear density at the resonance decay position r; po = 0.16 fm~ 3 is 
the normal nuclear density and a ~ 0.16 for the p and a ~ 0.12 for the uj meson [73] . 
The parametrization ( jTUl) may be employed also at much higher collision energies 
(e.g. FAIR and SPS) and one does not have to introduce a cut-off density in order 
to avoid negative pole masses. Note that (fTUj) is uniquely fixed by the 'customary' 
expression M^{p N ) w M (l — (%Pn/Po) in the low density regime. 

The spectral function of the vector meson V for the mass M at baryon density p^ 
is taken in the Breit-Wigner form: 



A(Mo)-C 2 M 2 Ty(M, p N ) 

A V (M, PN ) -C X .- (M2 _ M ^ W))2 + (MrMM;P7v))2 • (") 
The factor C\ is fixed by the normalization condition for arbitrary p^\ 



J A v (M,p N )dM = 1, (12) 

where Mu m = 2 GeV is chosen as an upper limit for the numerical integration. 
The lower limit for the vacuum spectral function corresponds to the two-pion decay, 
M m i n = 2m- K , whereas for the in- medium collisional broadening case M m i n = 2m e —>■ 
with m e denoting the electron mass. Mq is the pole mass of the vector meson 
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Fig. 2. The spectral functions for the p and uj meson in the case of the 'collisional broaden- 
ing' scenario (upper part) and the 'dropping mass + collisional broadening' scenario (lower 
part) for nuclear densities of 0,l,2,3,5xpo as employed in the transport calculations (see 
text for details). 

spectral function which is Mq(pjv = 0) = M in vacuum, however, may be shifted in 
the medium for the dropping mass scenario according to Eq. (11 01) . 

The resulting spectral functions for the p and uj meson are displayed in Fig. [2] for 
the case of 'collisional broadening' (upper part) as well as for the 'dropping mass 
+ collisional broadening' scenario (lower part) for densities of 0,1,2,3,5 xp - Note 
that in vacuum the hadronic widths vanish for the p below the two-pion mass and 
for the uj below the three-pion mass. With increasing nuclear density p^ elastic and 
inleastic interactions of the vector mesons shift strength to low invariant masses. 
In the 'collisional broadening' scenario we find a dominant enhancement of strength 
below the pole mass for the p meson while the uj meson spectral function is drastically 
enhanced in the low- and high-mass region with density (on expense of the pole- 
mass regime). In the 'dropping mass + collisional broadening' scenario both vector 
mesons dominantly show a shift of strength to low invariant masses with increasing 
Pn- Qualitatively similar pictures are obtained for the meson but quantitatively 
smaller effects are seen due to the lower effect of mass shifts and a substantially 
reduced (pN cross section which is a consequence of the ss substructure of the 
meson. Since the </> dynamics turn out to be of minor importance for the dilepton 




dropping mass 
+ collisional broadening 
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spectra to be discussed below we discard an explicit representation. 

The 'family' of spectral functions shown in Fig. 1 allows for a sufficient flexibility 
with respect to the possible scenarios outlined above. A comparison to dilepton data 
is expected to provide further constraints on the possible realizations. 

2.2 Off-shell propagation 

The underlying concept of all 'on-shell' semi-classical transport models (VUU, BUU, 
QMD, UrQMD etc.) is based on a Monte-Carlo solution of the transport equations 
for the phase-space densities of hadrons which are obtained in first order gradient 
expansion of the Wigner transformed Kadanoff-Baym equations using the quasipar- 
ticle approximation (QPA). The latter (QPA) assumes that the particles have narrow 
(i.e. 5-like) spectral functions which don't change their shape during the propaga- 
tion. This concept works sufficiently well for weakly interacting particles with a very 
long lifetime but becomes inadequate for short-lived (broad) resonance states of high 
collision rate. 

A first attempt to dynamically account for in-medium effects of vector mesons (i.e. 
a broadening of their spectral function) using the standard on-shell transport model 
showed a couple of severe problems (cf. the study by Effenberger et al. in Ref. [93]). 
The most serious one is that e.g. a p meson with a broad spectral function in the 
medium does not regain its vacuum spectral shape when propagating out of the 
medium The reason is that dynamical changes of the spectral function during 
propagation are only included by explicit collisions with other particles in conven- 
tional on-shell models (cf. the discussion in context with Fig. 2 below). 

In order to allow for a proper phase-space description of broad resonances with vary- 
ing spectral functions the off-shell transport dynamics (OSTD) has been developed in 
Refs. [&il65] . The main idea of the OSTD is to discard the quasiparticle approxima- 
tion in the first order gradient expansion of the Wigner transformed Kadanoff-Baym 
equations. This leads to generalized transport equations which contain an additional 
'backflow' term that vanishes in the on-shell QPA, however, survives in the off-shell 
case and becomes very important for the dynamics of broad resonances [SlKJo] (see 
below). A first application of the OSTD to dilepton production has been presented 
already in 2001 by one of the authors [48J. It has been shown that using OSTD the 
broad spectral function of vector mesons merges to the vacuum shape by propagation 
out of the dense medium. 

Here we briefly describe the implementation of OSTD in the HSD model: A vector 
meson V - created at some space-time point x at density p^ - is distributed in 
mass according to its spectral function A V (M, p N ) (and the available phase space). 
While propagating through the nuclear medium the total width of the vector meson 
Ty(M, pn) changes dynamical and its spectral function ffTTl) is also modified 
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according to the real part of the vector meson self energy ReY7 et as well as by the 
imaginary part of the self energy (T v ~ —ImY7 et /M) following 

2 M 2 T* 

A V (M, PN ) =C X -~ {M 2 _ M 2 _ ReE re t) 2 + {Mnr ( 13 ) 

which is the in-medium form for a scalar boson spectral function. The spectral func- 
tion has to change during propagation and to merge the vacuum spectral function 
when propagating out of the medium. 

In the OSTD the general off-shell equations of motion for test particles with momen- 
tum Pi, energy £j at position Xi read [5~4"f65] 



dX,, 



dt 1 — C(j) 2s i 



- - , e 2 - P 2 - Ml - ReEffi 
2Pi + V P , ReY$ + — 2-5 ^ V P , r , 



(0 



(0 



(14) 



dP 



dt 



1 — C(i) 2ej 



e 2 - P 2 



Ml 



ret 



r (0 



(15) 



dsj 
dt 



1 — 2£j 



P 1 



M 2 - ReE$ dY 



(0 



dt 



<9f 



(16) 



where the notation implies that the function is taken at the coordinates of the 
test particle, i.e. F (i) = F(t, X t (t), P^t), e^t)). In Eqs. (jHjl -fllEj) ReE ret denotes the 
real part of the retarded self energy while T = —ImY7 et stands for the (negative) 
imaginary part in short-hand notation. Apart from the propagation in the real po- 
tential ~ ReY7 et /2e the equations f|T4|) - f[T6]) include the dynamical changes due to 
the imaginary part of the self energy ImY7' et ~ — MT V with T v from ([6]). 

In Eqs. ( 1141) - (1161) a common factor (1 — C^)) -1 appears, which includes the energy 
derivatives of the retarded selfenergy, 



1 

271 



dReZ r $ 



+ 



e 2 - P 2 



Ml 



JfeE# 



(0 



(0 



(17) 



This common factor leads to a rescaling of the 'eigentime' of particle % but does not 
change the trajectories as demonstrated in Fig. 3 of Ref. [65]). According to the model 
studies in [65] this prefactor is negative for invariant masses M below the pole mass 
Mq and positive above which implies that the factor 1/(1 — Cj) is less than 1 for M < 
M and larger than 1 for M > M (as stated in [63]). This yields a time dilatation for 
masses M < Mo in their phase-space propagation. In order to examine the possible 
effect on dilepton spectra we have repeated the model calculations discussed in Fig. 
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3 of |65j including additionally the decay of each 'testparticle' according to the width 
r/(2Mj) in the eigen frame of the particle. Furthermore, Eqs. (22), (23) in Section 2.5 
of the present manuscript have been adopted to calculate the influence on dilepton 
spectra by integrating the differential decay rate in time. As initial condition 'broad 
vector quasiparticles' close to the center of a 'nucleus' have been assumed whereas 
the other parameters have been taken the same as in Fig. 3 of [65]. The resulting 
'dilepton spectra' show a slight enhancement with respect to the reference calculation 
for masses M < M when including the correction factor. This results from a slightly 
larger 'lifetime' of the particles for M < M - when including the rescaling of the 
'lifetime' by the correction factor - and leads to a longer dilepton radiation time. The 
small differences obtained can be traced back to the size of the correction factor for 
M < Mq for the model case in Ref. [65] which is less than 12%. Accordingly we find 
the 'enhancement' in the dilepton yield to be less than 12% for the model scenario, 
too. Since the explicit energy dependence of the retarded selfenergy is highly model 
dependent we will discard the terms (1171) in the following study, i.e. assume = 0. 

The interpretation of the equations of motion (|14p - ([TBI) becomes particularly trans- 
parent if T is independent on the 3-momentum P. Then, using M 2 = P 2 — ReY7 et 
as an independent variable instead of the energy P = e, Eq. (fl6|) turns to 



dM[ M 2 - M 2 dY {i 
dt f(j) dt 



(18) 



for the time evolution of the test-particle i in the invariant mass squared [64,65j. 
Now the deviation from the pole mass, i.e. AM 2 = M 2 — Mq, follows the equation 



d A w2 AM 2 d~ d ( AM 2 \\ . , 

—AM = -r <-> — In =0, 19 

dt r dt dt\\Y)J 

which expresses the fact that the off-shellness in mass scales with the total width. 

In order to demonstrate the importance of off-shell transport dynamics we present 
in Fig. [3] the time evolution of the mass distribution of p (upper part) and uj (lower 
part) mesons for central C+C collisions (b=l fm) at 2 A GeV for the dropping mass 
+ collisional broadening scenario (as an example). The l.h.s. of Fig. [3] corresponds 
to the calculations with on-shell propagation whereas the r.h.s. show the results for 
the off-shell dynamics. 

As seen from Fig. [3] the initial p and u mass distributions are quite broad even for a 
small system such as C+C where, however, the baryon density at 2 A GeV may reach 
(in some local cells) up to 2p . The number of vector mesons decreases with time due 
to their decays and the absorption by baryons (pN — > ttN or pN — > ttttN). Most of 
the p mesons decay /disappear already inside the 'fireball' for density pn > 0. Due 
to the 'fireball' expansion the baryon density drops quite fast, so some amount of p 
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Fig. 3. Time evolution of the mass distribution of p (upper part) and u (lower part) mesons 
for central C + C collisions (b=l fm) at 2 A GeV for the dropping mass + collisional broad- 
ening scenario. The l.h.s. of Fig. [3] correspond to the calculations with on-shell dynamics 
whereas the r.h.s. show the off-shell HSD results. 

mesons reach the very low density zone or even the 'vacuum'. Since for the off-shell 
case (r.h.s. of Fig. [3]) the p spectral function changes dynamically by propagation 
in the dense medium according to Eqs. (14)-(16), it regains the vacuum shape for 
p — > 0. This does not happen for the on-shell treatment (l.h.s. of Fig. [3]); the p 
spectral function does not change its shape by propagation but only by explicit 
collisions with other particles. Indeed, there is a number of p's which survive the 
decay or absorption and leave the 'fireball' with masses below 2m n . 

Accordingly, the on-shell treatment leads to the appearance of p mesons in the vac- 
uum with M < 2m n , which can not decay to two pions; thus they live practically 
'forever' since the probability to decay to other channels is very small. Indeed, such 
p's will continuously shine low mass dileptons which leads to an unphysical 'enhance- 
ment/divergence' of the dilepton yield at low masses (note, that the dilepton yield 
is additionally enhanced by a factor ~ 1/M 3 ). 
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The same statements are valid for the u mesons (cf. lower part of Fig. [3]): since the 
u} is a long living resonance, a larger amount of cj's survives with an in-medium like 
spectral function in the vacuum (in case of on-shell dynamics). Such uu's with M < 
3m T can decay only to 717 or electromagnetically (if M < m n ). Since such unphysical 
phenomena appear in on-shell transport descriptions (including an explicit vector- 
meson propagation) an off-shell treatment is mandatory. 

The off-shell equations of motion (fl4|) - ([16!) are the present standard for the HSD 
transport approach and have also been adopted by the Rossendorf group in Ref . [57] . 



2. 3 Vector meson production in the nuclear medium 



Apart from the explicit off-shell propagation the production of mesons is also modi- 
fied in the nuclear medium in line with their in-medium spectral function. We men- 
tion that in-medium mass shifts of particles have been traditionally incorporated in 
on-shell transport approaches by parametrizing the vacuum cross sections as a func- 
tion of the invariant energy yfs and the threshold y/so (cf. [2011211176] ). In- medium 
modifications of cross sections then in first order have been incorporated by a shift 
of the 'pole' threshold (which corresponds to the pole mass of the spectral function) 
V^o = y/so + — Mo due to a lack of microscopic calculations for the corresponding 
in-medium transition matrix elements. 

In order to account for the in-medium effects in production cross sections, we model 
the vector meson production cross sections in NN and ttN reactions in the following 
way: The total cross section ctnn^vnn(s, pn) (similar o~ n N->VN(s, pn)) is 



o~nn->vnn(s, pn) 



M„ 



dM 



da 



NN—fVNN 



[s, M, p N ) 



dM 



(20) 



The mass differential cross sections are approximated by 

da NN ^ V NN{s, M, p N ) 



dM 



o-° NN ^ VNN {s, M, p N ) ■ A(M, p N ) 



J A(M lPN )dM 
J A(M,p N )dM 



(21) 



where A(M,pn) denotes the meson spectral function (iTTj) for given total width Ty 
(JH1); M max = y/s — Itxin is the maximal kinematically allowed invariant mass of the 
vector meson V and is the mass of the nucleon. In Eq. (12~TT) cr NN _ fVNN (s, M, p^) 
is the parametrization (based on phase space) of the vacuum cross section [21] with 
a modified 'pole' threshold \/s^(pn) = 2niN + Mq(pn) instead of ^/sq = 2ijin + M Q 
in the vacuum. Note, that the physical threshold for the p production, e.g., in iVTV 
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Fig. 4. The production cross sections for the channels NN — » VX as a function of the invari- 
ant energy y/s for densities pj\r = and 3/3o- The solid lines represent the parametrizations 
of inclusive NN — > VX cross sections in free space while the dash-dotted lines stand for 
the 'collisional broadening' scenario; the dashed lines show the inclusive cross sections in 
the 'dropping mass + collisional broadening' scenario. The vertical thin blue line displays 
the threshold for meson production by string formation and decay (y/s th = 2.6 GeV) in 
case of baryon-baryon channels. The experimental data [77|7 8 88J are shown for exclusive 
pp — > Vpp (dots) and inclusive pp — ► VX (squares) vector meson production where X 
stands for two baryons and further mesons. 




2 3 456789 10 2 3 456789 10 



s [GeV] s w [GeV] 

Fig. 5. The production cross sections for the channels ttN — > VN (V = p, to) as a function 
of the invariant energy yfs for densities pn = and 3po- The solid lines represent the 
parametrizations of inclusive ttN — > VX cross sections in free space while the dash-dotted 
lines stand for the 'collisional broadening' scenario. The dashed lines show the inclusive 
cross sections in the 'dropping mass + collisional broadening' scenario. The experimental 
data |77|88j are shown for exclusive ttN — * VN (dots) and inclusive irN — * VX (squares) 
vector meson production. The vertical light blue line shows the threshold for meson pro- 
duction by string formation and decay (\/s th = 2.3 GeV) as implemented in HSD for 
meson-baryon channels. 



reactions is denned as \/s th = 2m N + M min , i.e. for the broadening scenario \/s th — > 
2m n- Thus, the formula fl2T|) is used to model the vector meson production in NN 
reactions in the vacuum and in the medium. 
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The results for the channels NN — > VX are displayed in Fig. H] as a function of the 
invariant energy yfs for densities p^ = and Pn = 3po- The solid lines represent 
the parametrizations for inclusive (NN — > VX) cross sections in free space while 
the dash-dotted lines stand for the 'collisional broadening' scenario; the dashed lines 
show the inclusive cross sections in the 'dropping mass + collisional broadening' 
scenario. The light vertical blue line stands for the threshold for meson production 
by string formation and decay (y/s th = 2.6 GeV) for baryon-baryon channels. The 
squares show the experimental data [771l78|88j for the inclusive NN — > VX vector 
meson production (here X stays for all final particles including two baryons and 
further mesons). The dots correspond to the experimental data for the exclusive 
pp — > Vpp reaction. 

The string formation and decay reasonable matches the inclusive (NN — > VX) data 
for p (squares) taken from Refs. [77|I78|I88| . It is seen that the excitation function 
of the vector meson production cross section in the 'collisional broadening' scenario 
is basically smeared out close to threshold. Only when incorporating additionally a 
dropping mass the thresholds are shifted down in energy such that the production 
cross sections become enhanced in the subthreshold regime with increasing nuclear 
density p N . 

The actual results for the channels nN — ► VN (V = p, to) are displayed in Fig. [5] as 
a function of the invariant energy \fs for densities p^ = and p^ = 3p - Here the 
solid lines represent the parametrizations of nN — > VX cross sections in free space 
while the dash-dotted lines stand for the 'collisional broadening' scenario. The dashed 
lines show the inclusive cross sections in the 'dropping mass + collisional broadening' 
scenario. The experimental data are displayed for exclusive tcN — > VN (dots) and 
inclusive irN — > VX (squares) vector meson production channels. The vertical light 
blue line stands for the threshold for meson production by string formation and decay 
(\fs th = 2.3 GeV) as implemented in HSD for meson-baryon channels. It is seen that 
the string decay reasonable matches the inclusive (ttN — > VX) data for p and to 
(squares) taken from Refs. [77f88] . 

One observes from Fig. [5] that in the 'collisional broadening' scenario the excitation 
function of the vector meson production cross section is only smeared out close to 
threshold due to the broader vector meson width at finite density. As in case of NN 
reactions the thresholds are shifted down in energy when incorporating addition- 
ally a dropping mass. Since the collisional broadening is low for the <fi meson the 
modifications of the production cross section for 0's stay very moderate (not shown 
explicitly). 

A note of caution has to be added here because the present parametrizations for 
vector meson production in the nuclear medium are not controlled by microscopic 
studies on the various production channels for the case of finite nuclear density p^. 
A consistent microscopic description of in-medium vector meson production cross 
sections thus is urgently needed. Furthermore, it should be noted that the experi- 
mental determination of a p meson production cross section is problematic since this 
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would imply to measure, e.g., the full invariant mass range of 7r + 7r~ pairs in a p- 
wave final state. This is a rather complicated task because the two-pion decay width 
T nn approaches zero close to the two-pion threshold where the data are dominated 
by pions in a relative s-wave. Accordingly experimental data have to be considered 
simultaneously with the operational definition for a 'p-meson'. 

We also mention that in the early study [15] we have considered the possible contribu- 
tion from the Dalitz decay of the iV(1520) resonance: iV(1520) — ► pN — > e + e~N. At 
those times this process has been treated perturbatively, i.e. without explicit propa- 
gation of the iV(1520). We obtained from our qualitative estimates that the dilepton 
radiation from the iV(1520) in the hot and dense medium is strongly reduced due 
to a 'melting' of the iV(1520). Since the width and pole position of the iV(1520) is 
still quite uncertain we do not include the contribution of the iV(1520) explicitly and 
rather rely on the findings in [15] that the iV(1520) decay contribution to dileptons 
should become small in a dense medium. 



2.4 'In-medium' extension of the LUND string model 

The particle production from baryon-baryon collisions above \/s th ~ 2.6 GeV and 
from meson-baryon collisions above yfs th — 2.3 GeV are modeled in HSD via string 
formation and decay. The actual realization is the LUND string-fragmentation model 
Fritiof [67] which is known to describe inelastic hadronic reactions in a wide energy 
regime at zero baryon density. In the standard version of Fritiof the resonances are 
produced according to the vacuum Breit-Wigner spectral function (by default non- 
relativistic) with a constant width. Moreover, the Breit-Wigner shape is truncated 
symmetrically around the pole mass, \M — Mq\ < 5, with 5 chosen 'properly' for each 
particle such that no problems are encountered in the particle decay chains [67J. 

However, in order to study in-medium effects at energies close to the individual 
thresholds \fs th one has to incorporate the relativistic in-medium spectral functions 
(as described in Section 2.1 for the vector mesons) also consistently in the string 
fragmentation model. Such an extension of the Fritiof model is performed in this 
study for the first time and the resonance production in the medium via Fritiof 
is treated with the full relativistic Breit-Wigner spectral functions including the 
density dependent self-energy and in-medium width (depending on mass and baryon 
density). Also the truncation of the spectral function in mass is removed, i.e. the 
resonance mass is chosen within the physical thresholds from M min to M max . As 
before the total energy and momentum conservation holds strictly in the extended 
Fritiof model. 

We mention that this extension of the Lund model - which allows to implement 
any shapes of spectral functions into the fragmentation scheme - has been used not 
only for vector mesons but also for the A resonance in order to incorporate the 
mass-dependent width T^(M^) (in line with the low energy part of HSD, where 
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the parametrization from Ref. [16] has been used for Ia(Ma)), which leads to a 
hardening of the high mass dilepton spectra from the A Dalitz decay. Note, that 'in- 
medium' effects - i.e. collisional broadening or/and dropping masses - are accounted 
only for the vector mesons (p, u>,(f>) in the present study. 



2. 5 Time integration ( or 'shining ') method for dilepton production 



Since dilepton production is a very rare process (e.g. the branching ratio for the 
vector meson decay is ~ 10~ 5 ), a perturbative method is used in order to increase 
statistics. In the HSD approach (in this study as well as in our earlier investigations 
PT|44f45ll48f49j ) we use the time integration (or 'shining') method first introduced 
by Li and Ko [71]. The main idea of this method is that dileptons can be emitted 
during the full lifetime of the resonance R before its strong decay into hadrons or 
absorption by the surrounding medium. For example, the p° decay (with invariant 
mass M) to e + e~ during the propagation through the medium from the production 
time t = up to the final ('death') time tp - which might correspond to an absorption 
by baryons or to reactions with other hadrons as well as the strong decay into two 
pions - is calculated as 

dM ^ 1 ] 7 (ftc) AM 1 ' 



in the mass bin AM and time step At (in fm/c). In (I22p 7 is the Lorentz factor 
of the p-meson with respect to the calculational frame. The electromagnetic decay 
width is defined as 

A/f * 4 

r^ e+e "(M) = C Pl J_ (23) 



where C p = Y p ^ e e (M )/M . Here M is the vacuum pole mass, Mq is the in- 
medium pole mass which is equal to the vacuum pole mass for the collisional broad- 
ening scenario, however, is shifted for the dropping mass scenario according to Eq. 

<m- 

The time integration method allows to account for the full in-medium dynamics of 
vector mesons from production ('birth') up to their 'death'. We note that by calcu- 
lating the dilepton emission only at the strong decay vertex (e.g. as in Refs. [oTIfoT] ) 
the dilepton rate (as well as the density dependence of the dilepton emission) is un- 
derestimated since a sizeable part of the emission history (e.g. before the absorption 
point by baryons or before the decay to pions) is lost. 
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2.6 e + e bremsstrahlung in elementary reactions 



The implementation of e + e~ bremsstrahlung from hadronic reactions in transport 
approaches has been limited so far to the dominant pn channel in the 'soft photon' 
approximation (with phase-space corrections) [2"Tll22f4Uf41ll42j or even been neglected 
completely [5Q155] - As we will demonstrate below this contribution should be reex- 
amined and revised accordingly. 

The soft-photon approximation (SPA) [79] is based on the assumption that the radi- 
ation from internal lines is negligible and the strong interaction vertex is on-shell. In 
this case the strong interaction part and the electromagnetic part can be separated, 
however, the cross section for dileptons has to be corrected [79] by reducing the 
phase-space for the colliding particles in their final-state. The phase-space corrected 
soft-photon cross section for the reaction 1 + 2 — ► 1 + 2 + e + + e~ can be written as 



da _ a 2 a(s) R 2 (s 2 ) ^ 



dyd 2 q T dM Qir 2 Mq 2 R 2 (s) 
R 2 (s) = \Jl - (mi + m 2 ) 2 /s, 

s 2 = s + M 2 - 2g 0V ^, 
s — (mi + m 2 ) 2 



o-(s 



2m\ 



where m\ is the mass of the charged accelerated particle, M is the dilepton invariant 
mass, qo the energy, qx the transverse momentum and y the rapidity of the dilepton 
pair. In (I24j) er(s) is the on-shell elastic cross section for the reaction 1 + 2 — ^ 1 + 2. 



In spite of the general limitation of the SPA - discussed in detail in Refs. [79|80] 
- Eq. f)24p has been widely used for the calculation of the bremsstrahlung dilepton 
spectra by different transport groups [2Tll22M40|l4"l~ir4 2] . Note, at those early times the 
applicability of the SPA for an estimate of dilepton radiation from NN collisions 
at 1-2 GeV bombarding energies has been supported by independent One-Boson- 
Exchange (OBE) model calculations by Schafer et al. [81] and later on by Shyam 
et al. [82] (cf. Fig. [6]). In these models the effective parameters have been adjusted 
to describe elastic NN scattering at intermediate energies. The models have then 
been applied to bremsstrahlung processes including the interference of different di- 
agrams for the dilepton emission from all charged hadrons. As shown in Ref. [81] 
the pp bremsstrahlung is much smaller than pn bremsstrahlung due to a destructive 
interference of amplitudes from the initial and final radiation. We note that 'gauge 
invariant' results have been obtained in Refs. [8~Tf8"2"] by 'gauging' the phenomeno log- 
ical form factors at the meson-baryon vertices. However, there are different schemes 
to introduce 'gauge invariance' in OBE models - as stressed by Kondratyuk and 
Scholten [M] - which lead to sizeably different cross sections. 



One has to point out that already in 1997 an independent study by de Jong et al. [83] 
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based on a full T-matrix approach - has indicated that the validity of the SPA for 
bremsstrahlung at intermediate energies of 1-2 GeV may be very questionable. 



e 1 e 



However, in Ref. [53] only the pp reaction has been considered (cf. r.h.s. of Fig. [H]); 
indeed, the bremsstrahlung in the full T-matrix approach is larger by a factor of 
about 3 than the corresponding SPA calculations (and OBE results). 
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Fig. 6. The e + e~ bremstrahlung frompn (l.h.s.) and pp (r.h.s.) channels. The dashed lines 
show the SPA results |79j . the dotted and dashed-dot-dotted lines correspond to the OBE 
calculations by Schafer et al. [H] and Shyam et al. [52] . respectively. The dash-dotted line 
(r.h.s.) displays the T-matrix result from de Jong et al. [83] while the red solid lines show 
the calculations by Kaptari et al. [69J. 

Recently, new covariant OBE calculations for dilepton bremsstrahlung have been 
performed by Kaptari and Kampfer [69] . The effective parameters for iVW scattering 
have been taken similar to [82J , however, the restoration of gauge invariance has been 
realized in a different way. As mentioned above, there are several prescriptions for 
restoring gauge invariance in effective theories including momentum-dependent form 
factors in interactions with charged hadrons [H] and the actual results very strongly 
depend on the prescription employed. The scheme in Ref. [69] is to include explicitly 
the vertex form factors into the Ward-Takahashi identity for the full meson-exchange 
propagators, which is different from the method used in Refs. [8Tf82] . 

The pn and pp bremsstrahlung differential spectra from Kaptari et al. [69] are shown 
in Fig. [6] by the solid red lines. One can see that the difference between the Kaptari 
model and the SPA and OBE [8TH82] results is up to a factor of 4 for pn reactions 
and about a factor of 3 for pp collisions at 1 GeV. Note, that the pp bremsstrahlung 
from [69] is on the same level as the T-matrix calculations from Ref. [83] (cf. Fig. 
El r.h.s.). Thus, according to the Kaptari model the contributions from pn and pp 
bremsstrahlung have been substantially underestimated in the previous transport 
model calculations [2T1[22] (accounting only for elastic scattering events with charged 
hadrons) . 



We stress that formula ff2Tj) is only applicable for the bremsstrahlung radiation from 
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elementary elastic scattering reactions (i.e. elastic iViV or ttN collisions). However, 
in the transport models of Refs. [56.57J it has been adopted also for inelastic NN 
reactions which faces the problem of double counting since the dilepton emission from 
inelastic NN reactions is dominantly taken into account by decays of the produced 
particles (i.e. A, rj, uj Dalitz decays, direct p,uj,<p decay etc.). On the other hand, 
electromagnetic radiation also emerges from inelastic channels since each charged 
particle (that is decelerated or accelerated in a collision) will radiate dileptons. Due 
to an interference of all amplitudes there is a partial cancellation between initial state 
amplitudes and final state amplitudes in case of equal charges due to a different sign 
in the acceleration. On the other hand two hadrons with opposite charges in the 
final state will also interfere destructively such that general estimates for a wide 
class of inelastic hadronic reactions cannot easily be performed. In this respect the 
strategy in Refs. [56.57J (to account for radiation from inelastic channels) is a very 
crude assumption whereas our approximation implies to disregard such additional 
radiation amplitudes in case of dominant electromagnetic decays which are taken into 
account explicitly. Due to a lack of microscopic calculations on the T-matrix level for 
radiation from channels with further charged particles in the final state we have to 
wait for respective experimental results to verify/falsify the different approximations. 



In the present study we will adopt the results from Ref. [69] for the elementary 
bremsstrahlung processes. As will be shown in the next Sections this will have a 
sizeable impact on our results for the dilepton radiation from AA collisions at 1 - 2 
A-GeV. 



We have to point out, furthermore, that in order to separate the bremsstrahlung 
(pp — ► ppe + e~) from a vector-dominance like dilepton production via the p-meson 
(pp — > ppp, p — > e + e~), we do not employ a vector-dominance formfactor when calcu- 
lating the bremsstrahlung. Thus, the dilepton radiation via the decay of the virtual 
photon (pp — > PP7*,7* — > e + e~) and the direct p decay to e + e~ are distinguished 
explicitly in the calculations. This is different from the calculations of the Tubingen 
group [17] which have no dynamical p meson but only a medium-dependent form- 
factor in the bremsstrahlung channel. 



In closing this Subsection we stress that new calculations of dilepton bremsstrahlung 
- based on a consistent full G— (or T— ) matrix approach - are urgently needed 
in order to clarify the unsatisfactory situation. From the experimental side high 
statistics measurements for multi-differential dilepton spectra from elementary pp 
and pn(pd) reactions are required to settle this open question for elastic as well as 
inelastic channels. 
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3 7r° and 77 production 



Since the low mass dilepton spectra are dominated by the Dalitz decays of it and 
77 an explicit control of their production cross section is mandatory for a proper 
interpretation of dilepton mass spectra. We recall that the early suggestion of the 
DLS Collaboration that the enhanced e + e~ production seen at the BEVELAC in 
C + C and Ca + Ca collisions might be due to an enhanced production of 77's was 
rejected by the TAPS Collaboration measuring the 77 cross section via the 7r°7 decay 
channel for the same systems [SS]. The TAPS argument has been confirmed by the 
transport calculations in Refs. [^|86] . 



In the last decade the experimental information of 77 production in pp and pn reac- 
tions - especially close to threshold - has been drastically improved such that more 
appropriate parametrizations of these cross sections have to be implemented in trans- 
port in order to stay in line with the experimental progress. The presently available 
experimental information from exclusive pp and pn reactions is displayed in Fig. [7J by 
the full symbols as a function of the invariant energy above threshold and indicates 
that the pn cross section is larger than the pp cross section for 77 production by up 
to a factor of 6. Only at sufficiently high energy the initial isospin no longer plays 
any role due to multiple associated pion production. The solid and dashed lines in 
Fig. [7J show the novel parametrizations/results from HSD for pp and pn reactions 
which smoothly join the inclusive pN — > nX production from string decay at high 
energy and are compatible with the low energy exclusive pN — > npN data from Refs. 

The number of 7r°'s and 77's from HSD for 12 C + 12 C at 1 A-GeV and 2 A-GeV are 
shown in Figs. [8] and [9] as a function of the impact parameter b (lower blue lines). 
The average number of 7r's and 77's is given by, 



_ 2nJ^N x (b)bdb 

~ ly x -> \bmax ~ 70 ) \ zo ) 

710 



max 



where b max denotes the maximum impact parameter considered in the calculation 
and x = (tt°, 77). The corresponding results for these quantities are displayed in Figs. 
E] and M by the red squares as a function of b = b max and demonstrate that the 
average number of 7r°'s and 77's from the TAPS Collaboration in inclusive 12 C + 12 C 
reactions (shaded areas) are reasonably reproduced by the transport calculations. 
Note that the inclusive reaction cross section is given by a inc i = ^b 2 max . 

In line with Figs. M and [9] we show in Fig. [10] the average tt° and 77 multiplicities for 
i2 C + i2 q from o g tQ 2 A-GeV from the TAPS Collaboration [Eq (full squares) in 
comparison to the corresponding HSD results (stars). Since the agreement is accept- 
able (within error bars) we may proceed with the actual dilepton studies. 
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Fig. 7. The rj production cross section from pp and pn reactions as a function of the invariant 
energy above threshold. The solid and dashed lines represent the corresponding inclusive 
(pN — > r]X) parametrizations/results from the HSD transport approach, whereas the full 
symbols show the exclusive (jpN — * rjpN) experimental data from Refs. [87 88 89J. 
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Fig. 8. The number of 7r°'s (l.h.s.) and 77's (r.h.s.) from HSD for 12 C+ 12 C at 1 A- GeV as a 
function of the impact parameter b (lower blue lines). The red squares denote the average 
number of 7r's and 77's according to Eq. (|25H for b max = b and show that the average number 
of 7r°'s and 77's from the TAPS Collaboration in inclusive 12 C+ 12 C reactions (shaded areas) 
is approximately reproduced for b max « 6 fm. 



4 Dilepton production in pp and pd reactions 



Before coming to explicit results for nucleus-nucleus collisions we consider dilepton 
production in more elementary reactions like pp and pd in the energy range of interest. 
The results for e + e~ production in pp reactions are shown in Fig. [TT] at bombarding 
energies of 1.04, 1.27, 1.61, 1.85, 2.09 and 4.88 GeV in comparison to the data from 
the BEVALAC [92] including the DLS acceptance filter and mass resolution [92J. 
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Fig. 9. The same as Fig. for 12 C + 12 C reactions at 2 A-GeV. 
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Fig. 10. The average 7r° and rj multiplicities for 12 C + 12 C from 0.8 to 2 A-GeV from the 
TAPS Collaboration [91] (full squares) in comparison to the corresponding HSD results 
(stars). 



The various contributions from the different channels (in the HSD calculations) are 
indicated additionally (see color coding in the figure). At all energies the 7T° Dalitz 
decay dominates for invariant masses up to the pion mass whereas the contributions 
from 0.2 to 0.6 GeV vary drastically with energy. At 1.09 and 1.27 GeV here the 
A Dalitz decay dominates which is superseeded by the r\ Dalitz decay already at 
1.61 GeV. The pp bremsstrahlung contribution also gives a sizeable yield to the total 
dilepton cross section at the lower energies. At 2.09 GeV we find that the calculations 
overestimate the yield from the p and u decays which has to be kept in mind when 
comparing to nucleus-nucleus data at 2 A-GeV (see below). At 4.88 GeV all channels 
result from string excitation and decay (except bremsstrahlung) which is sufficiently 
well in line with the BEVALAC data. We note that at all energies the contribution 
from the direct Dalitz decay of heavy resonances R — > e + e~ N is suppressed compare 
to A — > e + e~N as pointed out before in Ref. [94J. 
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Fig. 11. The differential cross section for e + e~ production in pp reactions at bombarding 
energies of 1.04, 1.27, 1.61, 1.85, 2.09 and 4.88 GeV in comparison to the data from the 
BEVALAC |92j including the DLS acceptance filter and mass resolution [92J. The various 
contributions from the different channels (in the HSD calculations) are indicated addition- 
ally (see color coding in the legend). 

The mass-differential dilepton cross section for pd reactions is compared in Fig. [12] to 
the BEVALAC data at the same bombarding energies including the DLS acceptance 
filter and mass resolution [92] . In this case the relative contribution from the various 
Dalitz decays and vector mesons decays stays about the same as for pp reactions 
except the contribution from bremsstrahlung which at the lowest energies is of the 
same order as the A Dalitz decay. Note that the p and u production is overestimated 
at 2.09 GeV as in the case oi pp reactions. Furthermore, coherent production of rj,p 
and oj on the deuteron has been discarded in the HSD calculations since the latter 
channels show up only very close to the individual production thresholds and become 
negligible in A + A reactions. 
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Fig. 12. The differential cross section for e + e~ production in pd reactions at bombarding 
energies of 1.04, 1.27, 1.61, 1.85, 2.09 and 4.88 GeV in comparison to the data from the 
BEVALAC |92j including the DLS acceptance filter and mass resolution [92J. The various 
contributions from the different channels (in the HSD calculations) are indicated addition- 
ally (see color coding in the legend). 



In view of the very limited DLS acceptance and mass resolution the overall descrip- 
tion of these 'elementary' reactions is acceptable such that more complex reactions 
can be addressed. 
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5 Dilepton production in nucleus-nucleus collisions 



A major aim of the HADES Collaboration was to clarify the 'DLS-puzzle', i.e. to 
verify/falsify the DLS data [39] from the experimental side. Accordingly, the same 
systems have been reinvestigated at 1 A- GeV with the HADES detector in order to 
clarify the issue. The situation, however, is not as easy since the DLS and HADES 
acceptances differ significantly [39^6Tj. However, recently the HADES Collaboration 
has performed a direct comparison with the DLS data for C+C at 1 A GeV [39] 
by filtering the HADES data with the DLS acceptance [63]. Both measurements 
were found to agree very well [63]! Thus, the 'DLS puzzle' has been solved from the 
experimental side; there is no obvious contradiction between the DLS and HADES 
data! The major goal for the transport models now is to verify the solution of the 
'DLS-puzzle' from the theoretical side. Note that transport calculations allow for a 
direct comparison between the DLS and HADES measurements by employing the 
different FILTER routines to the same set of calculated events. 



5.1 Comparison to the DLS data 

We start with a reinvestigation of the DLS data employing the novel cross sections for 
7] production, the modified A Dalitz decay contribution as well as the additional (and 
enhanced) bremsstrahlung channels as described above. The results of our transport 
calculation are displayed in Fig. [13] for 12 C + 12 C at 1.04 A-GeV in case of 'free' 
vector-meson spectral functions (upper part) and in case of the 'dropping mass + 
collisional broadening' scenario (lower part) employing the DLS acceptance filter and 
mass resolution. 

In fact, the situation has substantially improved compared to the early studies and 
the missing yield in the 'free' scenario is now reduced to a factor of about 1.5 in the 
mass region from 0.25 to 0.5 GeV. This is due to slightly higher contributions from 
A and rj Dalitz decays and a significantly larger yield from bremsstrahlung channels. 

To see the effects of these modifications quantitatively, we show in Table 1 a com- 
parison of our 'new' results with the 'old' ones from Refs. [4"1"|T4"5] - indicated as 
'new' /'old' - for rj and A Dalitz decays, NN Bremsstrahlung channels as well as for 
the sum of all contributions for C+C at 1 A GeV for invariant masses of 0.2, 0.35 
and 0.5 GeV. Both results correspond to the 'free' scenario. 

As seen from Fig. [131 the bremsstrahlung yield now is similar in shape and magnitude 
as the A Dalitz decay contribution and enhanced by factor up to 5 due to the 
novel NN bremsstrahlung cross section from [59] arid accounting for the contribution 
from pp bremsstrahlung additionally to pn. The contribution of the pion-nucleon 
bremsstrahlung is quite small in the C+C system due to the limited energy available 
in meson-baryon collisions and the moderate rescattering rate of pions. 
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Fig. 13. Results of the HSD transport calculation for the mass differential dilepton spectra 
in case of 12 C + 12 C at 1.04 A-GeV in comparison to the DLS data [39]. The upper part 
shows the case of 'free' vector-meson spectral functions while the lower part gives the 
result for the 'dropping mass + collisional broadening' scenario. In both scenarios the DLS 
acceptance filter and mass resolution have been incorporated. The different color lines 
display individual channels in the transport calculation (see legend). 

As in our previous studies the rj Dalitz decay is subdominant. We obtain only a 
slight enhancement (up to 30%) of the rj yield at 1 A GeV compared to the early 
study [33] since modifications of the rj production cross section in A^A" reactions (cf. 
Section 3) are relevant at higher energies. The enhancement of the high mass A 
Dalitz decay yield is related mainly to the modifications in the LUND string model 
by implementation of the mass dependent A width in the string fragmentation (cf. 
Section 2.4) which enhances the population of high mass A resonances. Additionally, 
we use now the parametrization from Ernst et al. [22] for the A — > Ne + e~ decay. 

Medium modifications for the vector mesons turn out to yield an enhancement in 
the region of 0.4 < M < 0.5 GeV but are very moderate due to the light system 
C + C. Though a slightly better description of the DLS data is achieved it would 
be premature to claim the presence of in-medium effects on the vector mesons from 
these data. 
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Table 1 

The mass differential dilepton yield 4fe [/xb/GeV] for 77 and A Dalitz decays, NN 
bremsstrahlung as well as the sum of all contributions for C + C at 1 A GeV. The left num- 
bers in each column stand for the present results whereas the right numbers show the early 
results from Refs. |44|45] - 'new'/'old'. The last column presents the DLS experimental 
data |39j. Both calculations refer to the 'free' scenario. 



M [GeV] 


V 


A 


Br ems. NN 


sum 


exp. data 


0.2 


1.75 / 1.3 


7.05 / 5.5 


5.0 / 1.88 


14.2 / 4.7 


15.5 ±2.6 


0.35 


1.15 / 0.9 


1.25 / 0.4 


1.6 / 0.35 


4.1 / 1.57 


6.1 ± 1.1 


0.5 


0.12 / 0.11 


0.16 / 0.03 


0.22 / 0.04 


0.54 / 0.28 


0.54 ±0.2 



Additional information is provided by the m Ca ± 40 Ca system at 1.04 A- GeV where 
the measured DLS spectra are reproduced in a qualitatively and quantitatively simi- 
lar manner by the transport calculations (cf. Fig. [T4"l) as in case of the lighter system. 
Again the combined bremsstrahlung channels provide a contribution in the same 
order as the A Dalitz decay; the rj Dalitz decay remains subdominant but the ttN 
bremsstrahlung increases compared to the C + C system due to more frequent pion 
rescattering on nucleons. In-medium effects for the vector mesons can be indentified 
in the transport calculations but are hard to see in the total dilepton spectra. 

Is the 'DLS puzzle' solved? This question can only be answered in a convincing 
manner by comparison with the recent HADES data. 

5.2 Comparison to HADES data 

The same multi-differential dilepton spectra - used for comparison with the DLS 
data - are now filtered by the HADES acceptance routines and smeared with the 
HADES mass resolution [ST]. A comparison of our calculations for the 12 C ± 12 C 
system at 1.0 A-GeV with the HADES data from Ref. [63] is presented in Fig. [15] 
and demonstrates that the agreement between data and transport calculations is 
reasonable for the HADES data, too. In this case the differential dilepton spectrum 
is divided by the average number of 7r°'s which in experiment is determined by half 
of the average number of (tt + + The spectra look slightly different in shape 
due to the much higher acceptance at low invariant mass where the tt° Dalitz decay 
practically exhausts the dilepton spectrum. As in case of the DLS data the 77 Dalitz 
decay turns out to be subdominant and the A Dalitz decay to be of similar magnitude 
as the combined bremsstrahlung contribution. Again the effect of in-medium vector- 
meson spectral functions is hard to see in the total spectra. 

A comparison to the HADES mass-differential dilepton spectra for 12 C ± 12 C at 2 
A-GeV [61J is presented in Fig. [TSJ for the 'free' (upper part) and the in-medium 
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Ca+Ca, 1.04 A GeV 
in-medium effects: CB+DM 
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Fig. 14. The mass differential dilepton spectra from HSD in case of Ca + Ca at 1.04 
A-GeV in comparison to the DLS data [39]. The upper part shows the case of 'free' vec- 
tor-meson spectral functions while the lower part gives the result for the 'dropping mass 
+ collisional broadening' scenario. In both scenarios the DLS acceptance filter and mass 
resolution have been incorporated. The different color lines display individual channels in 
the transport calculation (see legend). 

scenario (lower part). At the higher bombarding energy the rj Dalitz decay provides 
the dominant contribution in the mass region from 0.2 to 0.5 GeV followed by A 
Dalitz decays and the combined bremsstrahlung channels. The mass region around 
0.75 GeV is overestimated in the 'free' scenario, whereas including in-medium spec- 
tral functions for the vector mesons the description of the data is improved (very 
similar to Ref. |57j ) due a shifting of strength from the vector- meson mass regime to 
lower invariant mass. However, the in-medium effects for the light C + C system are 
only very moderate. 

We now turn back to the system l2 C + 12 C at 1.04 A-GeV and ask the question if 
the agreement in Fig. [TB] might be accidental? To address this issue we investigate 
the dilepton rapidity spectra in the center-of-mass system for different cuts on the 
invariant mass. The results of our transport calculations are shown in Fig. [17] where 
the mass integrated rapidity distributions are displayed on the upper panel, a cut 
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Fig. 15. Results of the HSD transport calculation for the mass differential dilepton spectra 
- divided by the average number of 7r°'s - in case of 12 C + 12 C at 1.0 A-GeV in comparison 
to the HADES data [63J. The upper part shows the case of 'free' vector- meson spectral 
functions while the lower part gives the result for the 'dropping mass + collisional broad- 
ening' scenario. In both scenarios the HADES acceptance filter and mass resolution have 
been incorporated. The different color lines display individual channels in the transport 
calculation (see legend). 

on low mass dileptons (M < 0.15 GeV) is performed in the middle part while the 
lower part shows the rapidity distribution for a cut M > 0.15 GeV. The individual 
colored lines display the contributions from the various channels in the HSD calcu- 
lations. Since the experimental analysis for the dilepton rapidity spectra (involving 
the same cuts) is still in progress we still have to wait for the data to verify /falsify 
our predictions. 

We recall that very similar investigations have been performed for 12 C + 12 C at 2.0 
A-GeV by the HADES Collaboration [62]. Figure [18] shows their mass integrated 
rapidity distributions for dileptons from 12 C + 12 C at 2.0 A-GeV (top panel) in com- 
parison to the HSD calculations for the in-medium scenario. A cut on low mass 
dileptons (M < 0.15 GeV) is performed in the upper part, the middle panel corre- 
sponds to the cut 0.15 GeV < M < 0.55 GeV while the lower part shows the dilepton 
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Fig. 16. The mass differential dilepton spectra from HSD in case of C-\- C at 2 A-GeV in 
comparison to the HADES data [61] . The upper part shows the case of 'free' vector-meson 
spectral functions while the lower part gives the result for the 'dropping mass + collisional 
broadening' scenario. In both scenarios the HADES acceptance filter and mass resolu- 
tion have been incorporated. The different color lines display individual channels in the 
transport calculation (see legend). 

rapidity distribution for M > 0.55 GeV. As mentioned above, the rj Dalitz decay con- 
tribution becomes dominant in the interval 0.15 < M < 0.55 GeV at 2.0 A-GeV, 
whereas the contribution of bremsstrahlung decreases. In total we find the HSD cal- 
culations to provide a good description of the experimental rapidity distributions for 
all invariant mass ranges considered. 

Apart from dilepton rapidity distributions (for different mass bins) also transverse 
momentum spectra of e + e~ pairs provide additional information. For comparison 
with future HADES data we show the results of our HSD calculations in Fig. [19] 
for the dilepton transverse momentum spectra in the mass range M < 0.15 GeV 
(upper part), 0.15 < M < 0.55 GeV (middle part) and M > 0.55 GeV (lower 
part) for l2 C + 12 C at 1.0 A-GeV. For the lowest invariant mass bin the transverse 
momentum spectrum is entirely dominated by the ir° Dalitz decays, however, followed 
by the bremsstrahlung contribution at low p?- In the interesting mass range 0.15 
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Fig. 17. HSD predictions for the mass integrated rapidity distributions for dileptons for 
12 C + 12 C at 1.0 A-GeV (top panel). A cut on low mass dileptons (M < 0.15 GeV) is 
performed in the middle part while the lower part shows the dilepton rapidity distribution 
for a cut M > 0.15 GeV. The individual colored lines display the contributions from the 
various channels in the HSD calculations (see color coding in the legend). 

< M < 0.55 GeV the bremsstrahlung and A Dalitz contributions dominate which 
might provide a handle to establish the strength of the bremsstrahlung channels 
independently from the mass differential spectra. 

Whereas the data analysis for the dilepton px spectra for l2 C + 12 C at 1.04 A-GeV is 
still in progress we may gain related information in comparison to HADES data at 2.0 
A- GeV. Accordingly, the dilepton transverse momentum spectra from the in- medium 
HSD calculations for l2 C + 12 C at 2.0 A-GeV are shown in Fig. [20] in comparison 
to the data from the HADES Collaboration [62]. The upper part corresponds to the 
mass range M < 0.15 GeV/c 2 , the middle part to 0.15 < M < 0.55 GeV/c 2 and the 
lower part to M > 0.55 GeV/c 2 . Indeed, we find again a quite reasonable agreement 
with the HADES data on the transverse momentum distributions for practically all 
mass intervals. A slight overestimation of the experimental data at lower pt is seen 
for the mass interval M > 0.55 GeV/c 2 (lower part of Fig. [201) . This is due to the 
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Fig. 18. The mass integrated rapidity distributions for dileptons from HSD for 12 C+ 12 C at 
2.0 A-GeV (top panel) in comparison to the HADES data [62J. A cut on low mass dileptons 
(M < 0.15 GeV) is performed in the upper part, the middle panel corresponds to the cut 
0.15 GeV < M < 0.55 GeV while the lower part shows the dilepton rapidity distribution 
for M > 0.55 GeV. The individual colored lines display the contributions from the various 
channels in the HSD calculations (see color coding in the legend). 

fact that the 'dropping mass + collisional broadening' scenario for the vector mesons 
leads to a shift of the px spectra to the low px region. This effect is stronger for heavy 
systems with larger baryon densities and might be used to distinguish the different 
in-medium scenarios from the experimental side. 



5.3 Predictions at SIS energies 



The DLS data on e + e~ production from elementary reactions like pp and pd colli- 
sions have been taken in a rather small range of phase space with limited statistics. 
This is sufficient for 'pilot experiments' but decisive measurements - allowing for 
a distinction between different model scenarios - are still missing. To this end the 
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Fig. 19. HSD predictions for the dilepton transverse momentum spectra for 12 C + 12 C at 
1.0 A-GeV. The upper part corresponds to the mass range M < 0.15 GeV, the middle part 
to the interval 0.15 < M < 0.55 GeV and the lower part to M > 0.55 GeV. The individual 
colored lines display the contributions from the various channels in the HSD calculations 
(see color coding in the legend). 

HADES Collaboration is in the process of performing experiments with pp and/or pd 
reactions at GSI with improved statistics, mass resolution and phase-space coverage. 
These new experiments will also shed further light on the size of the bremsstrahlung 
channels in competition with hadron Dalitz decays. In order to prove/disprove the 
interaction scenarios employed here in HSD we present predictions for the future 
measurements within the presently known HADES detector accept ancjjj. 



We start with elementary reaction channels: Figure [21] shows the differential cross 
section for e + e~ production in pp (upper part), pn (middle part) and pd (lower part) 
reactions at bombarding energies of 1.25 GeV (l.h.s.), 2.2 GeV (middle) and 3.5 
GeV (r.h.s.) including the actual HADES acceptance filter and mass resolution (that 



1 Any change in experimental cuts or in the data analysis in future will easily be taken 
into account since the HSD 'raw data' will be made available to the HADES Collaboration. 
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Fig. 20. The dilepton transverse momentum spectra from HSD for 12 C + 12 C at 2.0 A-GeV 
in comparison to the data from the HADES Collaboration [62J. The upper part corresponds 
to the mass range M < 0.15 GeV, the middle part to the mass range 0.15 GeV < M < 
0.55 GeV and the lower part to M > 0.55 GeV. The individual colored lines display the 
contributions from the various channels in the HSD calculations (see line coding in the 
legend). 

differs in the magnetic settings for the different systems and bombarding energy). 
We note that the dilepton spectra are normalized to the number of ir° taken from 
the HSD calculations and not to (tt + + 7r~)/2 as for the heavy-ion reactions. This 
is important to point out due to the strong isospin asymmetry in pp,pn and pd 
reactions. 

It is seen from Fig.lSTlthat for low invariant masses (< 0.15 GeV) the 7r° Dalitz decays 
dominates which can be controlled additionally by direct pion measurements. At 1.25 
GeV the A Dalitz decay (dash-dotted blue line) is by far the dominant channel in 
pp collisions whereas the NN bremsstrahlung channels (dash-dotted orange line) 
overtake for pn collisions at this energy. In pd reactions the A Dalitz decay and the 
bremsstrahlung channels contribute with roughly the same magnitude in the entire 
invariant mass range. The rj Dalitz decay slightly contributes in the pd case due to 
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Fig. 21. The differential cross section for e + e~ production in pp (upper part), pn (middle 
part) and pd (lower part) reactions at bombarding energies of 1.25 GeV (l.h.s.), 2.2 GeV 
(middle) and 3.5 GeV (r.h.s.) including the HADES acceptance filter and mass resolution 
(see legend for the different color coding of the individual channels) . 

the internal Fermi motion of the nucleons in the deuteron. Discarding the low mass 
7T° decays we find a dominance of rj and A Dalitz decays at 2.2 and 3.5 GeV for 
pp, pn and pd collisions up to invariant masses of 0.6 GeV. Bremsstrahlung channels 
here are lower by more than an order of magnitude and are hardly to access at 
this energy. On the other hand the 7r°, 77 and A Dalitz decays may be controlled 
by direct hadronic measurements. In the vector- meson mass range (~ 0.8 GeV) the 
dilepton yield is clearly dominated by the direct p and u decays; a separation of these 
channels is possible by measuring additionaly the 7r°7 decay channel of the uo meson 
as in Ref. [75]. Note that the contribution from pion rescattering on the spectator 
nucleon in case of pd reactions (dotted light green lines) is down by roughly two 
orders of magnitude and will be hard to identify experimentally. 

We continue with nucleus-nucleus collisions. Since the modifications of the p-meson 
spectral function only leave moderate traces in the differential dilepton spectra for 
12 C + 12 C (cf. Section 5.2) we provide predictions for Ar + KCl at 1.75 A- GeV, a 
system which is currently analyzed by the HADES Collaboration and might allow for 
a closer view on the p-meson properties in a denser medium. The HSD predictions for 
the mass differential dilepton spectra - divided by the average number of (7r + + 7r~)/2 
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Fig. 22. The HSD predictions for the mass differential dilepton spectra - divided by the 
average number of 7r°'s taken as (tt + + vr~)/2 - in case of Ar + KCl at 1.75 A-GeV. The 
upper part shows the case of 'free' vector-meson spectral functions while the lower part 
gives the result for the 'dropping mass + collisional broadening' scenario. In both scenarios 
the actual HADES acceptance filter and mass resolution have been incorporated (see legend 
for the different color coding of the individual channels) . 

- in this case are presented in Fig. [22J The upper part shows the case of 'free' vector- 
meson spectral functions while the lower part gives the result for the 'dropping mass 
+ collisional broadening' scenario. In both scenarios the actual HADES acceptance 
filter and mass resolution have been incorporated. Since the standard Dalitz decays 
are roughly the same in both scenarios let us concentrate on the direct vector meson 
decay channels, i.e. the dashed red lines for the p meson and the double-dot - dashed 
purple lines for the us meson. The effect of the finite mass resolution can clearly be 
seen in case of the u meson for the 'free' scenario. The peak appears at the pole 
mass and is dominantly smeared by the experimental mass resolution (except from 
the moderate vacuum decay width). 

In the 'dropping mass + collisional broadening' scenario the shift of the ui strength 
to lower invariant masses is clearly seen in the calculations but will be very hard 
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to access in experiment. In this respect it might be more favorable to look for the 
7r°7 decay in comparison to the e + e~ decay in order to observe experimentally in- 
medium effects for the u meson. In case of the p meson the effects of a dropping 
mass + collisional broadening become clearly visible for invariant masses above 0.3 
GeV resulting in a significant enhancement (by more than a factor of two) in the 
mass differential dilepton spectra. Actually this in-medium enhancement is slightly 
lower than predicted in the (more simplified) on-shell quasi-particle approach in the 
nineties [23112 If 44|T4"5] but should be measured experimentally in case of high statistics. 



6 Summary 

In this study a detailed analysis of dilepton production in C + C and Ca + Ca at 1-2 
A-GeV as well as in elementary pp and pd reactions at 1-5 GeV has been presented 
within the microscopic HSD transport approach that incorporates a full off-shell 
propagation of the vector mesons in line with Refs. [64165] . A comparison of dilepton 
spectra from 'elementary' pp and pd reactions with the DLS data in the energy range 
from 1 to 4.9 GeV specifies the isospin dependence of the various channels and shows 
the adequacy/deficiency of the parametrization or model assumptions in the HSD 
approach. 

Our present HSD model includes the following novel developments with respect to 
the previous studies a decade ago [2T|23ll44ir4"5] : 

• off-shell dynamics for vector mesons and explicit (nonperturbative) propagation 
of vector mesons also at SIS energies 

• improved cross sections for rj production in pp and pn reactions 

• extension of the LUND string model for the production of particles according to 
'in-medium' spectral functions 

• (enhanced) bremsstrahlung contributions according to Ref. [69] which proofs vital 
for a solution of the DLS puzzle. 

Furthermore, we have analyzed the tt° and rj production in C + C collisions at SIS 
energies. The HSD results for 7r° and rj productions have been found to be in a good 
agreement with the data from the TAPS Collaboration from 0.8 to 2 A- GeV. These 
studies pin down the dominant contributions from 7r° and 77 Dalitz decays in the 
low mass dilepton spectra and also to a large extent the contribution from A Dalitz 
decays which are closely linked to the pion multiplicity. 

On the vector meson side different scenarios of in-medium modifications of p's and 
cu's, such as collisional broadening and additionally dropping vector meson masses (cf. 
Fig. [2]), have been investigated and the possibilities for an experimental observation 
of in-medium effects in AA reactions has been discussed. We find that within the 
extended off-shell transport approach we achieve a simultaneous description of the 
differential dilepton spectra for 12 C+ l2 C at 1.04 A- GeV from the DLS Collaboration 
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as well as the HADES Collaboration which points towards a consistency of the 
independent measurements within error bars. In general the experimental spectra 
are better described by the in-medium scenarios for the vector mesons, however, the 
modifications of the spectra relative to the 'free' scenario are very moderate for the 
light C + C systems. 

In order to shed some light on perspectives for heavier systems we have provided 
detailed predictions for Ar + KCl at 1.75 A-GeV as well as elementary pp and pd 
reactions at 1.25 GeV, 2.2 GeV and 3.5 GeV. The latter systems are addressed also 
by the HADES Collaboration and in future will shed light on the contribution of 
bremsstrahlung channels in elastic and inelastic channels at about 1 A- GeV and 
the in-medium modifications of the p-meson at 1.75 A-GeV in the heavier nucleus- 
nucleus system. We stress that new calculations of dilepton bremsstrahlung - based 
on a consistent full G— (or T— ) matrix approach - are urgently needed from the 
theoretical side in order to clarify the presently open situation. 
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